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Abstract 

The //-capture process on 3 He leading to a neutron, a deuteron and a /i-ne- 
utrino in the final state is studied. Three-nucleon Faddeev wave functions 
for the initial 3 He bound and the final neutron-deuteron scattering states are 
calculated using the BonnB and Paris nucleon-nucleon potentials. The nuclear 
weak current operator is restricted to impulse approximation. Large effects 
on the decay rates of the final state interaction are found. The comparison to 
recent experimental data shows that the inclusion of final state interactions 
drastically improves the description of the data. 

PACS numbers: 
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I. INTRODUCTION. 



Recently considerable progress has been achieved in the calculations of three nucleon 
(3N) continuum states. It is now possible to generate exact 3N bound |IJ and scattering |2| 
states for any realistic nucleon - nucleon interaction even with inclusion of a three-nucleon 
force (3NF). This opens now the possibility to study the interaction of electromagnetic 
and/or weak probes with the 3N systems ( 3 He or 3 H nuclei) without introducing into the 
analysis the uncertainties due to inadequate approximate 3N states. Only using such exact 
states nuclear dynamics can be tested with such probes and important information can be 
gained on the corresponding hadronic current operators. 

Elastic and inelastic electron scattering on 3 He ( 3 H) as well as photodisintegration of 
3 He or pd capture are prominent examples of such processes and have been studied since 
many years with the hope to get insights into 3N bound state wave functions and into the 
hadronic current operator. As was shown in a series of recent papers on electromagnetic 
processes || || || a very important and unavoidable ingredient for their analysis is the 
exact treatment of the interaction among the three nucleons in the continuum states. It is 
also important, that the three-nucleon wave functions should be based on realistic nuclear 
forces. In this paper we would like to study the importance of final state interactions (FSI) 
between the nucleons in muon capture on 3 He. 

There are three final channels following the capture of a negative muon on 3 He: 



H~ + 3 He -> 3 H + u 
[T + 3 He -> d + n + (1) 
\l~ + 3 He — > p + n + n + v^. 



The capture rate for the H channel has been extensively studied using the elementary 
particle method and the impulse approximation |J J7|. When removing the nuclear uncer- 
tainties by using accurate three-body bound wave functions, this reaction offers the possi- 
bility to extract the induced pseudoscalar coupling constant with nearly the same precision 
as from the capture process by a free proton [?]]. 

The old calculations of the total decay rates for the dm/ M (two-body breakup of 3 He) and 
pnnz/^ (three-body breakup of 3 He) channels performed in |J using the impulse approxima- 
tion and 3N bound and scattering states generated within the Amado model with separable 
1 So and 3 Si two-nucleon interactions showed that scattering effects are large. The resulting 
total rates were in agreement with some of the limited and rather inaccurate experimental 
results while disagreed with others 0. 

Recently the first measurement of energy spectra for deuteron and proton leaving 3 He 
after nuclear muon capture leading to two-body and three-body breakup, respectively, has 
been reported [§] 0. The partial capture rates were compared in |J to simple plane wave 
impulse approximation calculations yielding fair agreement with the measured proton energy 
spectrum but underpredicting the measured rate of deuteron production by a large factor. 
It is the aim of the present paper to study if the inclusion of FSI by using an exact n- 
d scattering state can account for this discrepancy. In section II we describe our way to 
fully include FSI. In section III predictions for decay rates obtained with two realistic NN 
interactions: Bonn B [10] and Paris [[LlJ are shown and compared to the experimental values. 



We conclude in section IV. 
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II. THEORETICAL FORMALISM. 



Fig.l. depicts the kinematics of the reaction. Our basis of the muon capture on 3 He forms 
the more fundamental weak- interaction capture process on hydrogen: fi~ + p — ► n -\- v^. 
The initial state 

| % > — | \^ > | \l>3 He m P > 

is composed of the atomic K-shell muon wave function | ^ s M > with spin projection and 
the 3 He state | \l/3He m P > with spin projection m and four momentum P. We choose the 
lab system with P = 0. The transition leads to the final state 

| / > = \v fl s v > | # (_) P' > . 

where is the four momentum of the neutrino and s v its spin projection. The state 
| \{r(~) P > is the interacting neutron-deuteron state with overall four momentum P . 

The corresponding S-matrix element Sjj is given in first order perturbation theory and 
assuming the Fermi form for the interaction Lagrangian by 

S fi = i J d 4 x < f | C(x) \ i>= (2) 

= z(2^) 4 8(P f ~Pi)^* 

* < # H P' | I A (0)| ^3 He m P > * <v s u |4(0) | * s M > 

G 

V2 



z(2tt) 4 6{P f - Pi) -= L A I A . 



The leptonic matrix element L\ is known and is expressed in terms of the corresponding 
Dirac spinors for the neutrino u(u, s u ) and the muon u(jl, s^) by 

u(?,Sv)lx(l - 7bM/*,*m) . (3) 



(2nY 



In the nuclear matrix element neither the current operator I A (0) nor the nuclear wave func- 
tions are equally well under control. In the following we restrict ourselves to the impulse 
approximation (IA) and neglect the exchange current effects. The single nucleon current 
operator is parametrized in terms of weak-interaction formfactors gj' 



i A (o) = J2Y,J d P d p' flt (^ «, t Mp, s) * 



(2tt) s / 

* M - 2mgJ) 1 X + gXip + P') X + ftVf + ^V} * 

* r_u(p', s')a(p', s', t') , (4) 

where a, are standard nucleon creation and anihilation operators, r_ is the isospin lowering 
operator and k = p' — p. For the nucleon weak formfactors we use the standard values 
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9l^)= , 2 - , (5) 

v 1 ^ 0.71 

9 m = - 



2mgX{k 2 ) 
9^) 



0.71 [GeV] 2 < 

-1.262 

fc 2 ~ 

1.19 [GeV] 2 - 

-3.7 



0.71 [GeV] 5 

-2mgf(k 2 ) 
k 2 + m| 



with = 138.13 MeV. 

Nuclear wave functions are generated by the nonrelativistic Schrodinger equation with 
realistic nuclear forces. Therefore to be consistent, the nuclear-current should also be chosen 
nonrelativistically. After a nonrelativistic reduction of (4) and introducing standard Jacobi 
momenta p, q one gets for I A defined in Eq.(2) 

I A = < #(-) | i A | ^3 Hc m > , (6) 

where the momentum space matrix elements of i A are 

< M\ i A (Q) I P'q' > = 50 - $8(q> ~ (q~ |q))I A (?,Q) , 



*(q,Q) = T^yAaYi- - £-) - + g%2m)i(axi>) + 

(zty) 6 m 2m 2m 

with tx = |i? + q ; V — —Q. The treatment of the final state follows || Q. For the 
convenience of the reader we repeat the most important steps. 

The final scattering state | > = \^nJ > is decomposed into Faddeev components 
\Fi> 



,(-) 1 



l*# > = (1 + P) I Fx > , (8) 

where P is a sum of cyclical and anticyclical permutations of three nucleons. The Faddeev 
component obeys the equation 

\F 1 >= \® nd > + G^t^P \F 1 >, (9) 

and the state | $ nc j > describes the free relative motion of the final nucleon and the deuteron. 
Inserting (H) and (|9|) into (||) yields 

I M = IpwiAS + ^rescatt j (10) 
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where IpwiAS (symmetrized plane wave impulse approximation) corresponds to the case 
that in Eq.(6) < | is replaced by 4^ < $ nc j | (1 + P) and therefore no interactions 
between the outgoing nucleon and the deuteron are present. For the Irescatt term, which 
contains all rescattering, < | is replaced by < F\ \ PtiGo(l + P). The rescattering 



term can be written as 



r rescatt = -4<$ nd |p|^>, (ii) 



with 

| U» > = hGoil + P)i"(Q) | ^3 Hc > + hGoP | U» > . (12) 

This integral equation has the same integral kernel which one also finds in the 3N continuum 
and the same numerical methods can be used to solve it for any NN interaction. We solve 
Eq.fll2[) in a partial wave decomposition and in momentum space. For details we refer to 
[14J] |T5| The decay rate follows from S/, in a standard way |I7| and for the capture 



process with an unpolarized initial state and without polarization of the outgoing particles, 
it is given by 

,ll(2m'a) 3 f 8 G 2 



dT = (2tt) 2 --^ L / dPdp d dpJ{P -v-p d - Pn)(2n) 8 — J2 E I L aI A I (13) 

oll(2m'a) 3 /• _ ^^, „, 8 G 2 



m.c,, m n ,m.,i 



^ 2 ^ )2 LL[Zmal_ t ^^^^^ ^ £ | LaI a 



with P = (m M + m3He, 0) in the lab. system and m! = ^^"^ ■ The factors in (13) come 



from averaging over the spin of the initial particles, (22)' fr° m the states normalizations, 
((27r) 2 (27r) 8 ), and from the muon wave function (which is taken as an atomic K-shell wave 

1 

137.036 • 



function), ( ^ 2m ) , with a - — - 



III. RESULTS AND COMPARISON WITH EXPERIMENTAL DATA. 

To calculate the decay rates one has to evaluate the nuclear matrix elements I A . The 
initial 3N bound state of 3 He is always based on a 34 channel Faddeev calculation. In treating 
the 3N continuum we included all 3N partial waves with total NN angular momenta up to 
jmax = 2. In order to check the convergence of the results also calculations restricted to 
3N states with 1 S + 3 Si — 3 Di NN force components only and with j max = 1 have 
been performed. As NN force we use the Bonn B [l(| and Paris [jllj potentials. Let us first 
consider the case of a kinematically complete situation in the outgoing n-d channel. Since 
the outgoing channel contains 3 particles and the decay originates from an initial state of 
zero momentum one needs only 2 kinematical parameters to completely define the momenta 
of the outgoing particles. For these 2 parameters the neutrino and deuteron energies E u and 
Ed can be taken. The kinematically allowed range of neutrino energies is E u 6 [0, E^ ax ] with 

E max = (»y + rnsHc) 2 - (3m TO - \E D \) 2 
2{m ll + ms Re ) 
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For each neutrino energy there is a range of allowed deuteron energies as shown in Fig. 2. Each 
deuteron energy corresponds to a unique angle between neutrino and deuteron momenta. In 
Fig. 3 we show decay rates d ^ E for two values of the deuteron energy. There is a drastic 
increase in decay rates by a factor of ~ 200, when the interactions between the neutron 
and the deuteron in the final state are included. The increase is practically independent 
from the deuteron and neutrino energies. Fig.4 shows that the results obtained restricting 
to j m ax = 2 NN force components might be not yet fully converged. Nevertheless as is 
clear from Fig. 5 they are sufficient to analyse the existing data. The drastic effects of FSI 
show that any analysis of the experimental data should be based on realistic continuum 
wave functions. 

Up to now only one data set for the fi— capture on 3 He leading to the n + d + 
outgoing channel exists In order to compare our results to those data one has to 
numerically integrate the decay rate at every deuteron energy over the unobserved neutrino 
energy. This leads to the results presented in Fig. 5. It is cleary seen that while the PWIAS 
predictions drastically underestimate the data, the inclusion of FSI leads to a rather good 
agreement, which is fairly independent from the NN interaction used. Also the covergence 
in jmax is similar to the kinematically complete cases as can be seen in Fig. 6. 

IV. SUMMARY AND CONCLUSIONS. 

We calculated the decay rates for muon capture on 3 He leading to the n+d+z/^ channel. 
Realistic NN forces have been used (Bonn B and Paris potentials) and both the 3N bound 
state and the 3N continuum scattering wave functions were evaluated consistently solving the 
corresponding Faddeev equations. We used the most simple nonrelativistic single nucleon 
weak-current operator parametrized by the standard nucleon weak formfactors. Both in 
the kinematically complete und the uncomplete situations the predictions of the plane-wave 
impulse approximation for the final neutron-deuteron state drastically differ by a factor of 
about ~ 200 from the full result including final state interaction. The comparison to the 
only existing uncomplete data for this process reveals, that a good description of the data 
only results when the FSI is taken into account. The total failure of PWIAS in describing 
those data shows that any future analysis of this process should be performed with bound 
and scattering states generated consistently from a realistic 3N Hamiltonian. From the 
theoretical point of view mesonic exchange currents should also be added in future analysis. 
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Figure Captions 

Fig.l Kinematics for the \i~ + 3 He — > n + d + reaction. 
Fig. 2 Kinematically allowed neutrino and deuteron energies. 

Fig. 3 The decay rate dE ^ E at two values of the deuteron energy E d = 18 MeV and 
Ed = 29 MeV. The dashed and dotted lines are the PWIAS predictions for the Bonn 
B and Paris potentials, respectively. The solid and dashed-dotted lines are the Bonn 
B and Paris potential predictions when in the final state all NN force components up 
to j max = 1 are included. 

Fig. 4 The convergence in j max of the decay rate dE ^ dE at two deuteron energies as in Fig.3. 
The dotted, dashed and solid lines correspond to 1 5'o+ 3 5'i — 3 Di, j max = 1 and j max = 2 
Bonn B calculations, respectively. 

Fig. 5 The decay rate -^r for the case when neutrino is not observed. For the description 
of lines see Fig.3. The squares are experimental points from |I] ||. 

Fig.6 The convergence in j max of the decay rate For the decription of lines see Fig.4. 
The squares are experimental points as in Fig. 5. 
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FIGURES 
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FIG. 1. Kinematics for the fi~ + 3 He — ► n + d + v„ reaction. 
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FIG. 2. Kinematically allowed neutrino and deuteron energies. 
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FIG. 3. The decay rate d ^ d ^ E at two values of the deuteron energy = 18 MeV and 
Ed = 29 MeV. The dashed and dotted lines are the PWIAS predictions for the Bonn B and Paris 
potentials, respectively. The solid and dashed-dotted lines are the Bonn B and Paris potential 
predictions when in the final state all NN force components up to j max = 1 are included. 
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FIG. 4. The convergence in j max of the decay rate dEddE at two deuteron energies as in Fig. 3. 
The dotted, dashed and solid lines correspond to 1 Sq + 3 S± — 3 Di, j max = 1 and j max = 2 Bonn B 
calculations, respectively. 
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FIG. 6. The convergence in j max of the decay rate ^j^. For the decription of lines see Fig. 4. 
The squares are experimental points as in Fig. 5. 
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